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Research  proposal  1.  Approximation  schemes  for  stochastic  dif¬ 
ferential  equations  driven  by  a  fractional  Brownian  motion 

Foreword 

Consider  a  d-dimensional  stochastic  differential  equation  driven  by  an  m-dimensional 
fractional  Brownian  motion  BH  =  {B^,t  G  [0,  T]}  with  Hurst  parameter  H  G  (|,l)  of 
the  form 

pt  m  pt 

Xt  =  x0+  /  b(Xs)ds  +  J2  /  ^{Xa)dBl  (1) 

Jo  j=1  Jo 

The  stochastic  integrals  appearing  in  the  above  formula  are  path-wise  Riemann-Stieltjes 
integrals  constructed  using  Young’s  methodology  (see  [31]).  We  have  studied  numerical 
approximations  for  the  solution  to  equation  (1)  based  on  uniform  partitions  of  the  interval 
[0,  T\,  ti  —  i  —  0, . . . ,  n.  It  was  proved  by  Mishura  in  2008  (see  [26])  that  the  classical 
Euler  approximation,  denoted  by  Xn,  satisfies 

sup  \X\l  -Xi\  <Cenl~2H+e, 

0  <t<T 

for  any  real  number  e  >  0.  Moreover,  the  convergence  rate  n 1~2H  is  sharp  for  this 
scheme,  in  the  sense  that  n2n~ 1  [Y”  —  Xt\  converges  almost  surely  to  a  finite  and  non-zero 
limit.  This  shows  that  the  numerical  scheme  Xn  has  a  rate  of  convergence  different  from 
the  Euler-Maruyama  scheme  for  the  classical  Brownian  motion,  which  is  not  surprising 
because  for  H  —  \  the  sequence  converges  to  the  solution  of  the  corresponding  Ito 
equation,  and  for  H  >  |  we  are  dealing  with  Riemann-Stieltjes  pathwise  integrals. 

We  have  introduced  a  new  (modified  Euler)  approximation  scheme  for  equation  (1) 
that  takes  into  account  the  correlation  between  the  input  noise  dBf1  jdt  and  the  system 
Xt.  This  new  scheme  has  the  rate  of  convergence  y"1,  where 

if  \<H<  |, 

if  H  =  |,  (2) 

if  |  <  H  <  1 . 

In  particular,  the  rate  of  convergence  becomes  n~  2  when  H  is  formally  set  to  |,  which 
matches  the  rate  of  convergence  of  the  Euler  scheme  in  the  case  of  the  classical  Brownian 
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motion.  Furthermore,  we  have  derived  the  asymptotic  behavior  of  the  fluctuations  of 
the  error.  The  proof  of  these  results  is  based  on  the  techniques  of  Malliavin  calculus  or 
stochastic  calculus  of  variations. 

In  the  second  part  of  this  project,  we  have  studied  two  variations  of  the  time  discrete 
Taylor  schemes  for  rough  differential  equations  and  for  stochastic  differential  equations 
driven  by  fractional  Brownian  motions. 

Results 

The  results  obtained  in  this  research  project  are  included  in  following  papers: 

[1]  Y.  Hu,  Y.  Liu  and  D.  Nualart:  Rate  of  convergence  and  asymptotic  error  distribution  in 
Euler  approximation  schemes  for  fractional  diffusions,  Annals  of  Applied  Probability  26,  no. 
2,  1147-1207,  2016. 

Here  is  a  summary  of  the  main  results  proved  in  this  paper.  First  we  have  established 
the  following  result  on  the  rate  of  convergence  for  the  modified  Euler  approximation 
scheme,  that  we  denote  by  Xn : 

Theorem  1  Let  Xn  be  the  modified  Euler  approximation  scheme  for  Equation  (1).  We 
assume  that  b  :  Rd  — »  Rd  and  a  :  Rd  — »  Rdxd  are  three  and  four  times  continuously 
differentiable,  respectively,  with  bounded  partial  derivatives.  Then  for  any  p  >  1  there 
exists  a  constant  C  independent  of  n  (but  dependent  on  p)  such  that 

sup  E  [\Xf  —  Xt\p]p  <  C'yff1 . 

0  <t<T 

In  the  case  77  G  (|,  |]  we  have  derived  a  central  theorem  for  the  fluctuations  of  the 
error.  More  precisely,  the  process  7 n(X  —  Xn )  converges  in  law  to  the  solution  of  a  linear 
stochastic  differential  equation  driven  by  a  matrix-valued  Brownian  motion  independent 
of  B  as  n  tends  to  infinity.  A  formal  statement  is  included  in  the  next  theorem. 

Theorem  2  Let  77  €  (|,  |]  and  let  Xn  be  the  modified  Euler  approximation  scheme  for 
Equation  (1).  Assume  that  b  and  a  are  four  and  five  times  continuously  differentiable, 
respectively,  with  bounded  partial  derivatives.  Then  the  following  convergence  in  law  the 
space  C([0,  T])  holds  as  n  tends  to  infinity, 

{ln(Xt  —  Xf) ,  £  G  [0,  T]}  — >■  {Ut  ,t  G  [0,  T]}  ,  (3) 
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where  {Ut,t  G  [0, T]}  is  the  solution  of  the  linear  d-dimensional  stochastic  differential 
equation 


In  this  equation  W  is  a  matrix-valued  Brownian  motion  independent  of  the  fractional 
Brownian  motion  BH . 

The  proof  of  this  theorem  is  based  on  a  limit  theorem  for  weighted  sums  and  a  cen¬ 
tral  limit  theorem  for  quadratic  functionals  of  the  fractional  Brownian  motion,  which  is 
established  applying  the  fourth  moment  theorem  of  Nualart  and  Peccati  [29]. 

When  we  let  the  Hurst  parameter  H  converge  to  \  we  obtain  formally  the  classical 
results  for  the  Brownian  motion,  obtained,  for  instance  by  Kurtz  and  Protter  in  [22], 

In  case  H  G  (|,  1),  using  a  limit  theorem  in  Lp  for  weighted  sums,  we  have  derived 
the  IT-limit  of  the  normalized  error  n(Xt  —  Xff)  in  the  case  H  G  (|,  1). 

Theorem  3  Let  H  G  (|,  1)  and  a  and  b  are  four  and  five  times  continuously  differen¬ 
tiable,  respectively,  with  bounded  partial  derivatives.  Then 


n(Xt  -  X?)  — >•  Ut 


in  Lp{fV)  as  n  tends  to  infinity,  where  {Ut ,  t  G  [0,  T]}  is  the  solution  of  the  following  linear 
stochastic  differential  equation 


777.  777 

Ut=  Vb(Xs)Usds  +  J2  /  Vaj(Xs)UsdBi  +  J2  /  (VuV)PQd^ 

Jo  j=1  Jo  lJ=1  Jo 

XT~  A Vbb){Xs)ds +1-  r (Xba)(Xs)dBs  +  ^Yl  f(yajb){Xs)dBi , 

J  0  J  0  j ^  J  0 

where  Z  is  a  matrix-valued  generalized  Rosenblatt-type  process. 

We  have  also  derived  the  weak  rate  of  convergence  for  the  modified  Euler  scheme  Xn 
to  equation  (1). 

Theorem  4  Let  Xn  be  the  modified  Euler  approximation  scheme  for  equation  (1).  We 
assume  b  :  Rd  — >  Rd  and  a  :  Rd  — »  Rdxd  are  three  and  four  times  continuously  differen¬ 
tiable,  respectively,  with  bounded  partial  derivatives.  Then  for  any  function  f  :  — »  M, 
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which  is  three  times  continuously  differentiable  with  bounded  derivatives,  there  exists  a 
constant  C  independent  of  n  such  that 


sup 

o  <t<T 


E[/(A«)] 


E[/(A”)] 


<  Cn 


-1 


(4) 


If  we  assume  that  b  and  a  are  four  and  five  times  continuously  differentiable,  respectively, 
with  bounded  partial  derivatives,  then  for  each  t  E  [0,T], 


n|E[/(Xt)]  —  E[/(Xtn)]  j  converges  to 


a\T 


171  pi  pi  pi 


E 


-  ji=  i«/0  J 0  JO 


E<  Dt.Di 


Xf(Xt)AtTs(X^al)(Xs) 


u  —  s \2H  2\s  —  r\2H  2dudsdr 


r s(Xbb)(Xs)ds+  /  r s(Xba){Xs)dBs 
Jo 

m  ,.t 


E  /  r.{Voib)(x.)dBi 

j-i Jo 


as  n  — >  oc. 


[2]  Y.  Hu,  Y.  Liu  and  D.  Nualart:  Taylor  schemes  for  rough  differential  equations  and 
fractional  diffusions,  Discrete  and  Continuous  Dynamical  Systems  Series  B  21,  no.  9, 
3115-3162,  2016. 

In  this  paper,  we  study  two  variations  of  the  time  discrete  Taylor  schemes  for  rough 
differential  equations  and  for  stochastic  differential  equations  driven  by  fractional  Brow¬ 
nian  motions.  One  is  the  incomplete  Taylor  scheme  which  excludes  some  terms  of  an 
Taylor  scheme  in  its  recursive  computation  so  as  to  reduce  the  computation  time.  The 
other  one  is  to  add  some  deterministic  terms  to  an  incomplete  Taylor  scheme  to  improve 
the  mean  rate  of  convergence.  Almost  sure  rate  of  convergence  and  Lp-rate  of  convergence 
are  obtained  for  the  incomplete  Taylor  schemes.  Almost  sure  rate  is  expressed  in  terms 
of  the  Holder  exponents  of  the  driving  signals  and  the  ZT-rate  is  expressed  by  the  Hurst 
parameters.  Both  rates  involves  with  the  incomplete  Taylor  scheme  in  a  very  explicit  way 
and  then  provide  us  with  the  best  incomplete  schemes,  depending  on  that  one  needs  the 
almost  sure  convergence  or  one  needs  Lp-convergence.  As  in  the  smooth  case,  general 
Taylor  schemes  are  always  complicated  to  deal  with.  The  incomplete  Taylor  scheme  is 
even  more  sophisticated  to  analyze.  A  new  feature  of  our  approach  is  the  explicit  expres¬ 
sion  of  the  error  functions  which  will  be  easier  to  study.  Estimates  for  multiple  integrals 
and  formulas  for  the  iterated  vector  fields  are  obtained  to  analyze  the  error  functions  and 
then  to  obtain  the  rates  of  convergence. 
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Research  proposal  2.  Stochastic  partial  differential  equations 
driven  by  fractional  noises 

Foreword 

We  have  studied  the  linear  stochastic  heat  equation  (also  called  parabolic  Anderson  model) 
on  driven  by  a  general  multiplicative  centered  Gaussian  noise.  This  equation  can  be 
written  as 

dvL  1 

—  =  -A u  +  uW,  t  >  0,  x  G  Md,  (5) 

C/  L  Z 

with  initial  condition  u0(x),  where  u0  is  a  continuous  and  bounded  function.  In  the  above 
equation,  W  is  a  formal  Gaussian  field  that  has  a  covariance  of  the  form 


E 


Wt,XWS,y 


70  -  t)  A(x  -  y), 


(6) 


where  7  and  A  are  general  nonnegative  and  nonnegative  definite  (generalized)  functions 
satisfying  some  integrability  conditions.  The  product  appearing  in  the  above  equation 
(5)  can  be  interpreted  as  an  ordinary  product  of  the  solution  ut,x  times  the  noise  Wt>x 
(which  is  a  distribution).  In  this  case  the  evolution  form  of  the  equation  will  involve  a 
Stratonovich  integral  (or  path- wise  Young  integral).  The  product  can  also  be  interpreted 
as  a  Wick  product  and  in  this  case  the  solution  satisfies  an  evolution  equation  formulated 
using  the  Skorohod  integral.  We  have  considered  both  formulations. 

There  has  been  a  widespread  interest  in  the  model  (5)  in  the  recent  past,  with  several 
motivations  for  its  study.  First,  it  appears  naturally  in  homogenization  problems  for 
partial  differential  equations  driven  by  highly  oscillating  stationary  random  fields  (see 
[2,  19,  21],  On  the  other  hand,  equation  (5)  is  also  related  to  the  KPZ  growth  model 
through  the  Cole-Hopf’s  transformation.  In  this  context,  definitions  of  the  equation  by 
means  of  renormalization  and  rough  paths  techniques  have  been  recently  investigated 
in  [15,  18].  There  is  also  a  strong  connexion  between  equation  (5)  and  the  partition 
function  of  directed  and  undirected  continuum  polymers.  This  link  has  been  exploited  in 
[23,  30]  and  is  particularly  present  in  [1],  where  basic  properties  of  an  equation  of  type  (5) 
are  translated  into  corresponding  properties  of  the  polymer.  Finally,  the  multiplicative 
stochastic  heat  equation  exhibits  concentration  properties  of  its  energy.  This  interesting 
phenomenon  is  referred  to  as  intermittency  for  the  process  u  solution  to  (5)  (see  e.g 
[12,  13,  14]),  and  as  a  localization  property  for  the  polymer  measure  [8].  The  intermittency 
property  for  our  model  is  one  of  the  main  result  we  have  obtained. 

In  this  project  we  have  derived  existence-uniqueness  results,  Feynman-Kac  represen¬ 
tations,  chaos  expansions  and  intermittency  results  for  a  very  wide  class  of  Gaussian 
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noises  W  (including  in  particular  those  considered  in  [4,  10]),  for  both  Skorohod  and 
Stratonovich  type  equations  (5).  In  particular  we  have  obtained  some  lower  and  upper 
bounds  for  the  moment  of  oder  k,  for  all  k  >  2,  which  are  sharp  in  the  sense  that  they 
have  the  same  exponential  order  as  the  upper  bounds. 

A  second  part  of  the  project  has  been  devoted  to  study  the  one-dimensional  stochastic 
partial  differential  equation 

du  K,d 2u  .  ,T-  m 

di=  2d^  +  <J^W'  t~0,  XeR’  (7) 

where  by 

E[W(s,  x)W(t,  y)j  =  H(2H  -  1  )S0(t  -  s)\x  -  y\2H~ 2  (8) 

with  |  <  H  <  That  is,  W  is  a  standard  Brownian  motion  in  time  and  a  fractional 
Brownian  motion  with  Hurst  parameter  H  in  the  space  variable.  The  spatial  covariance  is 
not  a  locally  integrable  function  and  the  above  expression  is  formal.  A  rigorous  treatment 
of  this  covariance  requires  the  introduction  of  the  spectral  measure  /x(d£)  =  c\ 
where  C\,h  is  a  constant  depending  on  H.  The  standard  methodology  (see,  for  instance, 
[10])  to  handle  homogeneous  spatial  covariances  cannot  be  applied  here.  In  a  recent  paper, 
Balan,  Jobs  and  Quer-Sardanyons  [6]  proved  the  existence  of  a  unique  mild  solution  for 
equation  (7)  in  the  case  a(u )  =  au  +  b,  but  their  method  cannot  be  extended  to  general 
nonlinear  coefficients. 

We  have  been  able  to  establish  the  uniqueness  of  solutions  for  a  nonlinear  coefficient 
cr,  using  a  truncation  argument  inspired  by  the  work  of  Gyongy  and  Nualart  in  [17]  on 
the  stochastic  Burgers  equation  on  the  whole  real  line  driven  by  a  space-time  white  noise. 
For  the  existence,  we  have  applied  the  methodology  developed  in  the  work  of  Gyongy 
in  [16]  on  semi-linear  stochastic  partial  differential  equations,  which  consists  in  taking 
approximations  obtained  by  regularizing  the  noise  and  using  a  compactness  argument 
on  a  suitable  space  of  trajectories,  together  with  the  strong  uniqueness  result.  We  have 
also  established  the  Holder  continuity  of  the  solution  u  in  both  space  and  time  variables 
and  derived  upper  bounds  for  the  moments  of  the  solution  using  a  sharp  Burkholder’s 
inequality,  as  well  as  the  matching  lower  bounds  for  the  second  moment  by  means  of  a 
Sobolev  embedding  argument. 

A  third  part  of  the  project  has  been  devoted  to  study  the  position  of  the  high  peaks 
that  are  farthest  away  from  the  origin.  The  propagation  of  the  farthest  high  peaks  was 
first  considered  by  Conus  and  Khoshnevisan  in  [5]  for  a  one-dimensional  heat  equation 
driven  by  space-time  white  noise,  where  it  is  shown  that  there  are  intermittency  fronts 
that  move  linearly  with  time  as  at.  Namely,  for  any  fixed  p  G  [2,oo),  if  a  is  sufficiently 


6 


small,  then  the  quantity  supM>atE(|u(t,a;)|p)  grows  exponentially  fast  as  t  tends  to  oo; 
whereas  the  preceding  quantity  vanishes  exponentially  fast  if  a  is  sufficiently  large.  To  be 
more  precise,  the  authors  of  [5]  define  for  every  a  >  0, 

y(a)  :=  lim sup  -  sup  logM(|w(t, x)|p) ,  (9) 

t— 5>oo  t  \x\>at 

and  think  of  as  an  intermittency  lower  front  if  5? (a)  <  0  for  all  a  >  au-,  and  of  au 
as  an  intermittency  upper  front  if  S^{a)  >  0  whenever  a  <  au-  In  [5]  it  is  shown  that 
for  each  real  number  p  >  2,  0  <  au  <  <  oo,  and  when  p  —  2,  some  bounds  for  aL 

and  au  are  given.  In  a  later  work  by  Chen  and  Dalang  [3],  it  is  proved  that  when  p  —  2, 
there  exists  a  critical  number  a*  =  such  that  Jf{a)  <  0  when  a  >  a*  while  S^{a)  >  0 
when  a  <  a*  (this  property  was  first  conjectured  in  [5]).  On  the  other  hand,  using  a 
variational  approach  we  have  been  able  to  compute  the  exponential  growth  indices  for 
any  value  p  >  2. 

Inspired  by  the  aforementioned  works,  we  have  studied  the  multidimensional  stochastic 
heat  equation  (5)  driven  by  a  colored  noise,  both  in  space  and  time,  when  the  solution 
is  interpreted  in  the  Skorohod  sense.  Due  to  the  presence  of  the  time  covariance,  the 
propagation  speed  of  the  farthest  high  peaks  may  not  be  linear.  Thus,  in  contrast  to  (9), 
the  inequality  \x\  >  at  there  needs  to  be  replaced  by  \x\  >  at6t  for  some  suitable  function 
Qf  When  A  is  the  Riesz  kernel,  we  obtain  a  better  estimate  of  the  intermittency  lower 
front.  We  also  have  provided  explicit  formulas  for  S^{a),  for  arbitraty  p  >  2,  using  a 
variational  approach. 

Further  contributions  to  stochastic  partial  differential  equations  include  an  explicit 
formula  for  the  two-point  correlation  function  for  the  solutions  to  the  stochastic  heat 
equation  driven  by  a  space-time  white  noise,  the  analysis  of  fractional  in  time  equations 
and  the  spatial  asymptotics  for  the  stochastic  heat  equation  driven  by  a  Gaussian  noise 
which  is  whit  in  time  and  it  has  the  covariance  of  a  fractional  Brownian  motion  with 
Hurst  parameter  H  €  (|,  |). 

Results 

The  results  obtained  in  this  research  project  are  included  in  the  following  papers  whose 
main  results  are  described  below. 

[3]  J.  Huang,  Y.  Hu,  D.  Nualart  and  S.  Tindel:  Stochastic  heat  equations  with  general 
multiplicative  Gaussian  noises:  Holder  continuity  and  intermittency.  Electronic  Journal  of 
Probability  20  (2015)  1-50. 

Here  is  a  summary  of  the  main  results  proved  in  this  paper: 
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(i)  In  the  Skorohod  case,  the  mild  solution  to  equation  (5)  has  a  formal  Wiener  chaos 
expansion,  which  converges  in  L2(fl),  provided  7  is  locally  integrable  and  the  spectral 
measure  /i  of  the  spatial  covariance  satisfies  the  following  integrability  condition  (known 
as  Dalang’s  condition): 


Md0 
1  +  KI5 


<  00. 


(10) 


Moreover,  the  solution  is  unique.  The  proof  of  this  result  is  based  on  Fourier  analysis 
techniques,  and  covers  the  particular  examples  of  the  Riesz  kernel  and  the  Bessel  kernel 
considered  by  Balan  and  Tudor  in  [7] .  Our  results  also  encompass  the  case  of  the  fractional 
covariance  A(x)  =  nf=i  —  l)\xi\2Hi~2,  where  Hi>  \  and  condition  (10)  is  satisfied 

if  and  only  if  Yli=i  Hi  >  d  —  1.  This  particular  structure  has  been  examined  in  [20]. 


(ii)  Under  these  general  hypotheses  to  ensure  the  existence  and  uniqueness  of  the  solution 
of  Skorohod  type  one  cannot  expect  to  have  a  Feynman-Kac  formula  for  the  solution,  but 
one  can  establish  Feynman-Kac-type  formulas  for  the  moments  of  the  solution.  More 
precisely,  for  any  integer  k  >  2 


E  [ukt ;  J  =  E  B 


\~\uo(Bl  +  x)  exp 


i=  1 


y. 

^1  <i<j<k 


7 (s  —  r)A(Bls  —  B°r)dsdr  ) 

0  /  J 


(11) 


where  {FU;  j  =  1  is  a  family  of  d-dimensional  independent  standard  Brownian 

motions  independent  of  W.  The  formulas  we  have  obtained,  generalize  those  obtained  for 
the  Riesz  or  the  Bessel  kernels  in  [7,  20]. 

(iii)  Consider  the  following  more  restrictive  integrability  assumptions  on  7  and  /j:  There 
exists  a  constant  0  <  fd  <  1  such  that  for  any  i  G  1, 

0  <  7 (t)  <  Cp\t\~p 


and  the  measure  /i  satisfies 


f  Mg) 
L  1  +  kT-2" 


<  00. 


Under  these  assumptions,  we  have  derived  a  Feynman-Kac  formula  for  the  solution  u  to 
Equation  (5)  in  the  Stratonovich  sense.  An  immediate  application  of  the  Feynman-Kac 
formula  is  the  Holder  continuity  of  the  solution. 


(iv)  In  the  Stratonovich  case,  we  have  given  a  notion  of  solution  to  equation  (5)  using  two 
different  methodologies.  One  is  based  on  the  Stratonovich  integral  defined  as  the  limit  in 
probability  of  the  integrals  with  respect  to  a  regularization  of  the  noise,  and  another  one 


uses  a  path-wise  approach,  weighted  Besov  spaces  and  a  Young  integral  approach.  We 
show  that  the  two  notions  coincide  and  some  existence- uniqueness  results  which  are  (to 
the  best  of  our  knowledge)  the  first  link  between  pathwise  and  Malliavin  calculus  solutions 
to  equation  (5). 

(v)  Under  some  further  restrictions,  we  have  obtained  some  sharp  lower  bounds  for  the 
moments  of  the  solution.  Namely,  we  can  find  explicit  numbers  K\  and  n2  and  constants 
Cj,  Cj  for  j  =  1,  2  such  that 

C\  exp  (ci tKlkK2)  <  E  [\ut,x\k]  <  £2  exp  ( c2tKlkK 2) 

for  all  t  >  0,  1  6  Sd  and  k  >  2. 

[4]  Y.  Hu,  J.  Huang,  K.  Le,  D.  Nualart  and  S.  Tindel:  Stochastic  heat  equation  with  rough 
dependence  in  space.  Under  revision  for  the  Annals  of  Probability.  Arxiv:1505.04924vl.  57 
pages  long. 

Here  is  a  summary  of  the  main  results  proved  in  this  paper: 

(i)  Uniqueness  result  for  equation  (7):  Consider  the  space  Zf  formed  by  all  random  fields 
v(t,  x )  such  that 


IMU  :=  sup  ||u(t,  Olli^nxR)  +  sup  Nl_H  v{t)  <  00,  (12) 

te[o,T]  te[o,T]  2 

where  p  >  2  and 

•A/| -H,PV<J)  =  llu(^  •)  -  v(t>  ■  +  ^)llLP(nxK)l^|2'ff-2^  •  (13) 

Assume  that  0  is  differentiable,  its  derivative  is  Lipschitz,  cr(0)  =  0  and  for  p  >  the 

initial  condition  u0  is  nr  UP(M)  and 

[  \\u0(-)-uo(-  +  h)\\2LPm\h\2H-2dh<oo.  (14) 

Jr 

Then,  if  u  and  v  are  two  solutions  of  (7)  and  u,v  G  Zf,  for  every  t  G  [0,T]  and  x  G  M, 
u(t,  x)  =  v(t,  x),  a.s. 

(ii)  Existence  result  for  equation  (7):  Define 

J\fpu0(x)  =  f  I  |u0(x  +  h)  —  u0(x)|2|/i|_1_2/3(i/i 
\J  R 
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and  assume  that  for  some  f3  >  \  —  H  and  some  p  >  max(4H6_1,  p+H-1/2^  the  initial 
condition  u0  is  in  LP(M)  fl  L°°(M)  and 

supj\fpu0(x)  +  (  /  ||uo(0  -  uq(-  +  h)\\2LPm\h\2H~2dh]  <  oo  .  (15) 

xCM  \J  R  J 


Suppose  also  that  cr  is  differentiable  and  the  derivative  of  a  is  Lipschitz  and  <r(0)  =  0. 

Then  there  exists  a  solution  u  to  (7)  in  the  space  DXf  H'P,  where  H'P  is  the  family 
of  random  helds  u  satisfying 


sup  \\u(t,x)\\Lp(n)  +  sup 
te[o,T],xeR  te[o,T],xeM 


| u(t,  x  +  y)  —  u(t,  x ) 


l|p(n)l 


\y\2H  2dy  )  <  oo. 


(ii)  Moment  estimates  for  equation  (7):  Let  \  <  H  <  and  consider  the  solution  u  to 
equation  (7)  with  cr(u)  =  u  wit  initial  condition  is  uq{x)  =  1.  Let  n  >  2  be  an  integer, 
x  G  M  and  t  >  0.  Then  there  exist  some  positive  constants  Ci,C2,C3  independent  of  n,  t 
and  k  with  0  <  c\  <  C2  satisfying 

exp(cin1+^K1-sf)  <  E [un(t,x)\  <  c3exp  (c2n1+^K1~^t)  . 


[5]  Y.  Hu,  J.  Huang  and  D.  Nualart:  On  the  intermittency  front  of  stochastic  heat  equation 
driven  by  colored  noises.  Electronic  Communications  in  Probability  21,  no  21,  1-13,  2016. 

Here  is  a  summary  of  the  main  results  proved  in  this  paper.  Suppose  that  the  spectral 
measure  of  the  spatial  covariance  /j  satisfies  Dalang’s  condition.  For  any  real  number 
N  >  0,  we  define 

c„=  f  and  Dn  =  ^{(:\(\<N]  .  (16) 

J\€\>n  Kr 

On  the  other  hand,  we  assume  that  7  is  locally  integrable,  we  set  f*'y(s)ds  =  Tt.  The 
next  theorem  provides  an  upper  bound  for  the  upper  intermittency  front: 


Theorem  5  Let  u(t,x )  be  the  solution  to  equation  (5)  driven  by  a  noise  W  with  covari¬ 
ance  structure  (6).  Assume  that  uq  is  non-negative  and  supported  in  the  ball  Bm  =  {i£ 

:  |x|  <  M}.  Set  9t  =  \J D]ytC where 


Nt  =  inf  <  Y  >  0  :  CN  < 


(2tt) 


32(p  —  l)A2Tt 


(17) 
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Then,  for  any  integer  p  >  2  ,  we  have 

p(p)  :  =  inf  <  g  >  0  :  limsup  sup  logE[up(f,a;)]  <  0  >  <  1 .  (18) 

j  t-¥  OO  tOt  \x\>gt9t  \ 

We  have  derived  also  a  lower  bound  for  the  lower  intermittency  front,  but  only  in  the 
particular  case  where  the  spatial  covariance  is  the  Riesz  kernel. 

[6]  J.  Huang,  K.  Le  and  D.  Nualart:  Large  time  asymptotics  for  the  parabolic  Anderson  model 
driven  by  spatially  correlated  noise.  To  appear  in  Annals  of  the  Institut  Henri  Poincare. 
Arxiv:1509.00897v3  39  pages  long. 

In  this  paper  we  study  the  linear  stochastic  heat  equation,  also  known  as  parabolic 
Anderson  model,  in  multidimension  driven  by  a  Gaussian  noise  which  is  white  in  time 
and  it  has  a  correlated  spatial  covariance.  Examples  of  such  covariance  include  the  Riesz 
kernel  in  any  dimension  and  the  covariance  of  the  fractional  Brownian  motion  with  Hurst 
parameter  H  G  (|,  |]  in  dimension  one.  First  we  establish  the  existence  of  a  unique 
mild  solution  and  we  derive  a  Feynman-Kac  formula  for  its  moments  using  a  family  of 
independent  Brownian  bridges  and  assuming  a  general  integrability  condition  on  the  initial 
data.  In  the  second  part  of  the  paper  we  compute  Lyapunov  exponents,  lower  and  upper 
exponential  growth  indices  in  terms  of  a  variational  quantity.  The  last  part  of  the  paper 
is  devoted  to  study  the  phase  transition  property  of  the  Anderson  model. 

[7]  Y.  Huang,  K.  Le  and  D.  Nualart:  Large  time  asymptotics  for  the  parabolic  Ander¬ 
son  model  driven  by  space  and  time  correlated  noise.  Submitted  for  publication  Arxiv:16- 
7.00682vl  23  pages  long. 

In  this  paper  we  study  the  linear  stochastic  heat  equation  on  M  ,  driven  by  a  Gaussian 
noise  which  is  colored  in  time  and  space.  The  spatial  covariance  satisfies  general  assump¬ 
tions  and  includes  examples  such  as  the  Riesz  kernel  in  any  dimension  and  the  covariance 
of  the  fractional  Brownian  motion  with  Hurst  parameter  H  e  (|,  |]  in  dimension  one. 
First  we  establish  the  existence  of  a  unique  mild  solution  and  we  derive  a  Feynman-Kac 
formula  for  its  moments  using  a  family  of  independent  Brownian  bridges  and  assuming 
a  general  integrability  condition  on  the  initial  data.  In  the  second  part  of  the  paper  we 
compute  Lyapunov  exponents  and  lower  and  upper  exponential  growth  indices  in  terms 
of  a  variational  quantity. 

[8]  X.  Chen,  Y.  Hu,  S.  Tindel  and  D.  Nualart:  Spatial  asymptotics  for  the  parabolic  Anderson 
model  driven  by  a  Gaussian  rough  noise.  Submitted  for  publication.  Arxiv:1607.04092vl.  41 

pages  long. 
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The  aim  of  this  paper  is  to  establish  the  almost  sure  asymptotic  behavior  as  the  space 
variable  becomes  large,  for  the  solution  to  the  one  spatial  dimensional  stochastic  heat 
equation  driven  by  a  Gaussian  noise  which  is  white  in  time  and  which  has  the  covariance 
structure  of  a  fractional  Brownian  motion  with  Hurst  parameter  H  G  (|,  |)  in  the  space 
variable. 


[9]  L.  Chen,  Y.  Hu,  K.  Kalbasi  and  D.  Nualart:  Intermittency  for  the  stochastic  heat 
equation  driven  by  fractional  noise  in  time  with  H  G  (0,1/2).  Submitted  for  publication. 
Arxiv:1602.05617vl.  22  pages  long. 

This  paper  studies  the  one-dimensional  stochastic  heat  equation  driven  by  a  Gaussian 
noise  which  is,  with  respect  to  time,  a  fractional  Brownian  motion  with  Hnrst  parameter 
H  G  (0,1/2).  We  establish  the  Feynman-Kac  representation  of  the  solution  and  obtain 
both  lower  and  upper  bounds  for  the  LP(Q)  moments. 


[10]  L.  Chen,  Y.  Hu  and  D.  Nnalart:  Nonlinear  stochastic  time-fractional  slow  and  fast 
diffusion  equations  on  Md.  Submitted  for  publication.  Arxiv:1509.07763vl.  43  pages  long. 

This  paper  studies  the  nonlinear  stochastic  partial  differential  equation  of  fractional 
orders  both  in  space  and  time  variables: 


+  ^(— A)Q/2j  u(t,x)  =  //  p(u(t,x))W(t,x)  ,  t  >  0,  x  G 


where  W  is  the  space-time  white  noise,  a  G  (0,2],  /i  G  (0,2),  7  >  0  and  v  >  0.  Fun¬ 
damental  solutions  and  their  properties,  in  particular  the  nonnegativity,  are  derived  and 
proved.  The  existence  and  uniqueness  of  solution  together  with  the  moment  bounds  of  the 
solution  are  obtained  under  Dalang’s  condition:  d  <  2a  +  |  min(2q  —1,0).  In  some  cases, 
the  initial  data  can  be  measures.  When  (3  G  (0, 1],  we  prove  the  sample  path  regularity 
of  the  solution. 


[11]  L.  Chen,  Y.  Hu  and  D.  Nnalart:  Two-point  correlation  function  and  Feynman-Kac  for¬ 
mula  for  the  stochastic  heat  equation.  To  appear  in  Potential  Analysis.  Arxiv:1607.00682vl 

23  pages  long. 

In  this  paper,  we  obtain  an  explicit  formula  for  the  two-point  correlation  function  for 
the  solutions  to  the  stochastic  heat  equation  on  M.  The  bounds  for  p-th  moments  proved 
by  Chen  and  Dalang  in  [9]  are  simplified.  We  validate  the  Feynman-Kac  formula  for  the 
p-point  correlation  function  of  the  solutions  to  this  equation  with  measure-valued  initial 
data,  using  techniques  from  Malliavin  calculus. 
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Research  proposal  3.  Estimation  of  parameters  for  stochastic 
differential  equations  driven  by  a  fractional  Brownian  motion 

Foreword 

We  have  studied  fractional  Volterra  processes  X%  =  {Xi(t),t  >  0},  i  —  1, . , . ,  k,  of  the 
form 

Xi(t )  =  [  xjt  —  s)dBH(s),  t  >  0,  (19) 

Jo 

where  BH  is  a  fractional  Brownian  motion  with  Hurst  parameter  H  >  |,  and  Xi  :  [0,  oo)  — > 
M  are  measurable  functions  satisfying  suitable  integrability  conditions. 

The  special  case  of  A:  —  1  and  x± (u)  =  ae~9u,  with  a,  6  >  0,  corresponds  to  the 
fractional  Ornstein-Uhlenbeck  process.  More  generally,  our  model  includes  the  solution 
to  A:-th  order  stochastic  differential  equation  driven  by  BH ,  known  as  fractional  CARA: 
model,  which  generalizes  the  model  introduced  in  [24]: 

f  X^{t)  =  £ JlJ  OjXVKt)  +  aBH(t),  t  >  0 
\  X(0)  =  ...  =X(fc-1)(0)  =  0  ’  1  ’ 

where  the  superscript  ^  denotes  j-fold  differentiation  with  respect  to  t. 

Motivated  by  the  parameter  estimation  in  the  model  (20),  we  have  derived  general 
central  limit  theorems  for  functionals  of  the  process  (19).  The  results  are  applied  to  the 
problem  of  parameter  estimation  in  the  fractional  CARA:  model,  and  the  limit  matrix 
covariance  is  computed  explicitly  in  the  particular  case  k  =  2  and  under  some  restrictions 
on  the  parameters  90  and  0\ . 

In  a  second  part  of  this  research  direction,  we  have  derived  a  central  limit  theorem 
for  functionals  of  a  large  class  of  Gaussian  self-similar  processes.  We  have  also  studied 
self-similar  Gaussian  process  that  arise  form  stochastic  partial  differential  equations  with 
additive  noise.  For  these  processes  we  have  established  a  decomposition  in  law  and  a 
central  limit  theorem  for  the  Hermite  variations  of  the  increments.  Finally  we  have 
derived  a  general  Ito  formula  in  law  for  weak  symmetric  integrals  with  respect  to  the 
fractional  Brownian  motion,  we  have  studied  stochastic  differential  equations  with  power 
type  nonlinearities  and  we  have  derived  properties  of  the  derivative  of  the  self-intersection 
local  time  of  the  fractional  Brownian  motion. 

Results 

The  results  obtained  in  this  research  project  are  included  in  the  following  papers  whose 
main  results  are  described  below. 
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[12]  I.  Nourdin,  D.  Nualart  and  R.  Zintout:  Multivariate  central  limit  theorems  for  aver¬ 
ages  of  fractional  Volterra  processes  and  applications  to  parameter  estimation.  Accepted  in 
Statistical  Inference  for  Stochastic  Processes  19  no.  2,  219-234,  2016. 

Here  is  a  summary  of  the  main  results  proved  in  this  paper:  Consider  the  random 
vector  Ut  =  {I  i:r,  ■  ■  ■  >  E4,t),  where 


Ui,T  = 


(21) 


In  this  definition,  X,.  i  =  1, . . . ,  k,  are  the  fractional  Volterra  processes  defined  by  (19), 
fi,  i  —  1, . . .  ,  k,  are  real  mesurable  functions  satisfying 

f  fi(x)e~x2^2dx  =  0  and  f  f2(x)e~x2^2dx  <  oo,  (22) 

J  M  J  ]R 

and  aft)  =  y/ E[Xft)2].  The  second  condition  in  (22)  ensures  that  fi  can  be  expanded 
in  Hermite  polynomials,  namely 


/<  =  £«Mir,  with  ESo  l-ah  < 00 >  (23) 

1=0 

whereas  from  the  first  one  we  deduce  that  o  =  0. 

Then,  following  the  approach  developed  in  Nourdin,  Peccati  and  Podolskij  [28]  (see 
also  [27,  Chapter  7])  we  have  proved  the  following  result. 


Theorem  6  Let  qt  denote  the  Hermite  rank  of  fi,  that  is,  the  smallest  value  of  l  such 
that  the  coefficient  aiti  of  Hi  in  (23)  is  different  from  zero.  Set  g*  =  mui!^*.  g*  and 
assume  that  g*  >  2.  Consider  Ut  =  (Uitr,  ■  ■  ■  iU^t),  where  U^t  is  given  by  (21).  If 
H  e  (|,  1  —  A_)  and  if  the  functions  Xi  defining  Xj  satisfy  both 


\  9iV«?j 

xfu)xj(v)\  | v  —  u  —  a\2H~2dudv  ]  da  <  oo 


(24) 


and 


r)i  i  H(2H  —  1)  f  Xi(u)xfv)\v  —  u\2H~2dudv  e  (0,  oo), 
y  J[o,oo)2 


for  all  i,j  —  !,•••,  k,  then 


(25) 


UT  "  JV*(0,  A)  as  T  ->  oo, 


(26) 
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where  A  =  (A ij)i<i,j<k  is  given  by 


A ij  ^  ^  etijCLjjll 


H\2H  -l)1 


l=q-iVqj 


riWj 


(27) 


X 


'  [0,oo)2 


Xi(u)xj(v)\v  —  u  —  a\ltL  zdudv\  da. 


In  Theorem  6  we  must  divide  by  a  quantity  depending  on  t  in  (21),  namely  which 
is  not  very  convenient  for  applications.  In  this  sense  we  have  considered  the  random  vector 
VT  =  (Vi'T,  ■  ■  ■ ,  Vk.T)  given  by 


where  ry  is  given  in  (25),  and  we  have  proved  the  following  theorem. 

Theorem  7  Suppose  that  j\  =  P i  =  1, . . . ,  k,  are  real  polynomials  and  denote  by 
qi  the  Hermite  rank  of  Pi.  Set  q *  =  min \<i<kqi  and  assume  that  g*  >  2.  Consider 
Vt  =  iVi,Ti  ■  ■  ■ ,  Vk,T)  given  by  (28).  If  H  e  (|,  1  —  ^-)  and  if  the  functions  Xi  defining  Xt 
satisfy  (24),  (25)  as  well  as 

/  \xi(u)xi(v)\  ((u  A  v)  V  l)  |i>  —  u\2H~2dudv  <  oo,  (29) 

J  [0,oo)2 

then 

VT  ^  Nk( 0,  A)  asT  ->  oo,  (30) 

with  A  still  given  by  (27). 

[13]  D.  Harnett  and  D.  Nualart:  Central  limit  theorem  for  functionals  of  a  generalized  self¬ 
similar  Gaussian  process.  Arxiv:1508.02756vl.  Submitted  for  publication.  16  pages  long. 

Here  is  a  summary  of  the  main  results  proved  in  this  paper.  Let  A"  =  {Xtlt  >  0} 
denote  a  centered  self-similar  Gaussian  process  with  self- similarity  parameter  (5  e  (0, 1). 
Consider  the  following  conditions  on  the  function  (j)  defined  <f>{x)  =  E[XiXx]  for  x  >  1, 
where  a  G  (0,2/3]: 

(H.l)  (f  has  the  form  <f>{x)  =  —X(x  —  l)a  +  if(x),  where  A  >  0,  if{x)  is  twice-differentiable 
on  an  open  set  containing  [l,oo),  and  there  is  a  constant  C  >  0  such  that  for  any 
x  G  (1,  oo) 
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1.  \if'(x)\  A  Cx0-1; 

2.  \if)"(x)\  <  Cx^fx  —  l)a_1;  and 

3.  -0(1)  =  /h0'(l),  when  a  >  1. 


(H.2)  There  are  constants  C  >  0  and  1  <  v  <  2  such  that  for  all  x  >  2, 


1.  |(//(V)|  < 


l)-17 

f)a-2 


2.  W\x)\  < 


l)a-S 


if  a  <  1 
if  a  >  1, 

if  a  <  1 
if  a  >  1 


Under  these  conditions,  and  using  the  Fourth  Moment  Theorem  we  proved  the  follow¬ 
ing  result. 


Theorem  8  Suppose  a  self-similar  Gaussian  process  X  which  satisfies  ( H.l )  and  (H.2) 
above.  For  T  >  0  and  integers  n  >  2,  consider  the  sequence 

\nt\  — 1 

Hn  =  ~7=  f(Yj,n)i 

V  3=0 


_  A  Xj 


where  f  e  L2(M,  7)  has  Hermite  rank  d  >  2  and  7  =  7V(0, 1)  and  Yjjn  =  p- 

A  Xj  =  X(j+1yn  —  Xj/n.  Then,  if  a  <  %  —  the  sequence  {Fn,n  >  1}  converges  in 
distribution  to  a  Gaussian  random  variable,  with  mean  zero  and  variance  given  by  a2  = 
E7=yq<  where 


a2  =  2 ~qq\T  ^  (| m  +  l\a  +  \m  -  1|Q  -  2\m\a)q  . 

m£Z, 


[14]  D.  Harnett  and  D.  Nualart:  Decomposition  and  limit  theorems  for  a  class  of  self-similar 
Gaussian  processes.  Arxiv:1508.06641vl.  Submitted  for  publication.  16  pages  long. 

Here  is  a  summary  of  the  main  results  proved  in  this  paper:  Consider  a  centered 
Gaussian  process  {Xt,t  >  0}  with  covariance 


R(s,t)=E[XaXt] 


E 


Zt_rdB 5 


Z  q _ rdB(f 


(31) 


where 
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(i)  Bh  =  {Bfi ,  t  >  0}  is  a  fractional  Brownian  motion  with  Hnrst  parameter  H  E  (0, 1). 

(ii)  Z  =  {Zt,t  >  0}  is  a  zero-mean  Gaussian  process,  independent  of  BH ,  with  covari¬ 
ance 

E  [ZsZt\  =  ( s  +  t)~\  (32) 

where  0  <  7  <  2  H. 

In  other  words,  X  is  a  Gaussian  process  with  the  same  covariance  as  the  process  {  f{)  Zt_rdB[f ,  t  > 
0},  which  is  not  Gaussian.  Our  first  result  is  the  following  decomposition  in  law  of  the 
process  X  as  the  sum  of  a  fractional  Brownian  motion  with  Hurst  parameter  |  =  H  —  7 
plus  a  process  with  regular  trajectories. 

a 

Theorem  9  The  process  X  has  the  same  law  as  {y/fiB/  +  Yt,t  >  0},  where 


1  f°°  E7”1 

rfy)  Jo  i+^ 


dz , 


(33) 


B^  is  a  fractional  Brownian  motion  with  Hurst  parameter  a/2,  and  Y  (up  to  a  constant) 
is  the  process  introduced  by  Lei  and  Nualart  in  [25],  that  is,  Y  is  a  centered  Gaussian 
process  with  covariance  given  by 


where 


roc 

E  [YtYs]  =  X !  /  y~a~\  1  -  e~yt){l  -  e~ys)  dy, 

Jo 


Ai  = 


An 


roo  1-2  H 


r(^)T(2H  +  1)  sin(nH)  J0  1  +  r 


dr). 


Theorem  10  Let  q  >  2  be  an  integer  and  fix  a  real  T  >  0.  Suppose  that  a  <  2  —  K  For 
t  E  [0,  T],  define, 


Fn(t)  =  n 


\nt\  —  1 

E 

3= 0 


H„ 


(  AXj 


\l|AAj|| 


L2(Q) 


where  Hq(x )  denotes  the  qth  Hermite  polynomial  and  AXj  =  X(J+\yn  —  Xj/n.  Then 
as  n  — >  00,  the  stochastic  process  {Fn(t),t  E  [0,T]}  converges  in  law  in  the  Skorohod 
space  _D([0,  T]),  to  a  scaled  Brownian  motion  {oBt,t  E  [0,  T]},  where  {Bt,t  E  [0,  T]}  is  a 
standard  Brownian  motion  and  a  =  \fo ^  is  given  by 


o 


2 


Yq  dm  +  1I“  “  2\m\a  +  \m  ~  MaY  • 

mez 
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[15]  G.  Binotto,  I.  Nourdin  and  D.  Nualart:  Weak  symmetric  integrals  with  respect  to  the 
fractional  Brownian  motion.  Arxiv:1606.04046vl.  Submitted  for  publication.  21  pages  long. 

The  aim  of  this  paper  is  to  establish  the  weak  convergence,  in  the  topology  of  the 
Skorohod  space,  of  the  //-symmetric  Riemann  sums  for  functionals  of  the  fractional  Brow¬ 
nian  motion  when  the  Hurst  parameter  takes  the  critical  value  H  =  (4 1  +  2)_1,  where 
£  =  £(io)  >  1  is  the  largest  natural  number  satisfying  f()l  a2jis(da )  =  [  for  all  j  = 

0,  1.  As  a  consequence,  we  derive  a  change-of-variable  formula  in  distribution, 

where  the  correction  term  is  a  stochastic  integral  with  respect  to  a  Brownian  motion  that 
is  independent  of  the  fractional  Brownian  motion. 

[16]  J.  A.  Leon,  D.  Nualart  and  S.  Tindel:  Young  differential  equations  with  power  type 
nonlinearities.  Arxiv:1606.02258vl.  Submitted  for  publication.  25  pages  long. 

In  this  paper  we  give  several  methods  to  construct  nontrivial  solutions  to  the  equation 
dyt  =  c(yt)  dxt,  where  a;  is  a  7-Holder  Md-va]ued  signal  with  7  G  (1/2, 1)  and  o  is  a  function 
behaving  like  a  power  function  |£|K,  with  k  G  (0,1).  In  this  situation,  classical  Young 
integration  techniques  allow  to  get  existence  and  uniqueness  results  whenever  7(^+1)  >  1, 
while  we  focus  on  cases  where  7(rc  +  1)  <  1.  Our  analysis  then  relies  on  Zahle’s  extension 
of  Young’s  integral  allowing  to  cover  the  situation  at  hand. 

[17]  A.  Jaramillo  and  D.  Nualart:  Asymptotic  properties  of  the  derivative  of  self-intersection 
local  time  of  fractional  Brownian  motion.  To  appear  in  Stochastic  Processes  and  Their 
Applications.  Arxiv:1512.07219vl.  34  pages  long. 

Let  {Bt}t> 0  be  a  fractional  Brownian  motion  with  Hurst  parameter  |  <  H  <  1.  We 
prove  that  the  approximation  of  the  derivative  of  self-intersection  local  time,  defined  as 

ot£  =  (  [  p'e{Bt  -  Bs)As<At , 

Jo  Jo 

where  pe(x)  is  the  heat  kernel,  satisfies  a  central  limit  theorem  when  renormalized  by 

3  1 

£2~h.  We  prove  as  well  that  for  q  >  2,  the  y-tli  chaotic  component  of  a£  converges 
in  L 2  when  |  <  H  <  |,  and  satisfies  a  central  limit  theorem  when  renormalized  by  a 
multiplicative  factor  in  the  case  |  <  H  <  ^pEf. 
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